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Abstract 

The one solitary boson exchange potential (OSBEP) is used to evaluate 
observables of NN elastic scattering below pion threshold. In this approach, 
we use a nonlinear model of self-interacting mesons as a substitution for the 
commonly used phenomenological form factors. iViVdata support an empirical 
scaling law between the pion and other meson fields, which suggests a link 
to QCD and significantly reduces the number of parameters in the boson 
exchange potential. The analysis of np and pp observables distinguishes the 
model by its fit and few adjustable parameters. An outlook to apply OSBEP 
in vrA^ systems is given. 

PACS number(s): 13.75.Cs, ll.lO.Lm, 13.75. Gx, 21.30.-x 
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I. INTRODUCTION 



The notion of interacting elementary particles for low- and medium-energy nuclear 
physics is associated with definitions of potential operators which, inserted into a Lippmann- 
Schwinger equation, yield the scattering phase shifts and observables. In principle, this 
potential carries the rich QCD substructure consisting of quarks and gluons and thus may 
be deduced from some microscopic model. There are a number of models which explicitly 
refer to QCD and have gained remarkable success describing qualitative features of hadronic 
interactions Unfortunately, so far none of these models is able to reach the accuracy of 
phenomenological boson exchange or inversion potentials These models, however, do 

not contain any explicit reference to QCD and in the case of boson exchange models use ef- 
fective baryon and meson fields with phenomenological masses, coupling constants and form 
factors. It remains astonishing that, with these assumptions, they are able to account for a 
highly quantitative description of NN data below pion production threshold and thus have 
established themselves as the standard models to be used in nuclear physics. Furthermore, 
inversion and boson exchange models work equally well for meson-nucleon |^,^ and meson- 
meson interactions. This implies that the potentials remain valid at relative distances 
of ~0.3fm, which is much smaller than the rms radii of mesons and nucleons themselves 
and smaller than the QCD bag sizes. It is beyond any doubt that nucleons and mesons are 
genuine QCD objects and we expect their effects to become distinguishable within relative 
distances of ~1.5fm. In this context, it is common belief that phenomenological form factors 
effectively describe the actual QCD dynamics at short distances. 

To perform a step towards QCD inspired models, we attempt to replace the conventional 
form factors by a nonlinear meson dynamics using the one solitary boson exchange potential 
(OSBEP), which was developed recently by the Hamburg group 0. From the success 
of the empirical boson exchange potentials, it seems obvious that chiral symmetry is not 
dominant in A^A^ scattering below 300 MeV ||^ . None the less, also a phenomenological low- 
energy model should be inspired by concepts which ensure chiral symmetry conservation. 
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In this sense, we adopt structures from the hnear a model and develop a dynamics of self- 
interacting mesons. At this stage, we cannot circumvent chiral symmetry breaking by taking 
the nonlinearities, masses, and coupling constants as free parameters. Most important, the 
self-interaction is taken into account persistently at all instances. This is achieved by using 
meson fields which are quasiclassical analytic solutions of nonlinear field equations. Defining 
free meson operators, the quantization of these fields is done a posteriori. Finally, this model 
is utilized in the framework of a one boson exchange potential (OBEP), which closely follows 
the Bonn-B potential 

The benefit of this approach is the inclusion of nonlinear effects, leading to meson propa- 
gators of finite self-energy, which permits us to replace the form factors in conventional boson 
exchange potentials. Furthermore, an empirical scaling law was discovered which relates the 
pion mass and its self-interaction coupling constant with the self-interaction parameters of 
any of the other mesons used. Confirming our conjecture about reminiscent effects of the 
microscopic substructure subsumed in the empirical form factors, we interpret this as a hint 
for an underlying symmetry. An obvious benefit of the scaling law is the practical bisec- 
tion of the number of adjustable parameters. This is a different approach than pursued by 
the Bonn-CD |]TD[ or Reidlike Nijmegen potentials []TT|, which achieve perfect fits with an 



inflating number of parameters. 

A description of the theoretical framework, including technical details, can be found 
in 1^. Therein, we restricted the analysis to fit np SM95 phase shifts only. In the present 
work, we extend the potential to describe np as well as pp scattering and calculate scattering 
observables to be compared with the latest database compiled in SAID |jT^. Additionally, 
we show np and pp phase shift comparisons for Bonn-B 0, Nijm93 [|lT|, Paris |13|, OSBEP, 
and the analysis SM97 of Arndt et al. [p 



This paper is organized as follows. In Sec.|l|, we give the salient features of OSBEP. 
The fit of the model parameters to phase shifts is discussed in Sec.^ and thereafter, an 
extensive survey of np and pp scattering observables is given in Sec.^. An outlook for 
application of the nonlinear model to vrA^ scattering, together with a summary, is contained 



in Sec. 0. 



II. SOLITARY MESONS 

It is a common feature of chiral invariant models that spontaneous symmetry breaking 
leads to nonlinear terms in the mesonic part of a Lagrangian which can be interpreted as 



a self-interaction [|l5l . Because of this, a meson Lagrangian with the same structure as the 
linear cr model for all mesons in the OBE potential is assumed. Altogether, we consider the 
six mesons /5 = vr, r], p, u, a, 6, and a Lagrangian 

For mesons with nonzero spin the operator is a vector in Minkowski space. The parameter 
•p assumes 1/2 or 1 to distinguish odd and even powered nonlinearities and £int contains 
desirable couplings to nucleon and other meson fields. In chiral symmetric models, the 
self-interaction coupling constants Af and and the various meson masses are related by 
symmetry relations. This sounds intriguing but is not practical. In view of the ambiguities 
contained in Eq. ( |2.1| ) and, in particular in £int, it appears wise to restrict oneself first to 
a quantitative model which allows chiral symmetry breaking. In actual calculations, this 



implies the permission of free parameters in Eq. (|2.1D which are the coupling constants. 



physical masses, and the nonlinearities Af . After fitting the parameters to observables, we 
rely on their nature to effectively restore chiral symmetry a posteriori . 

A. Meson propagation 



The Lagrangian ( p. I] ) contains self-interacting mesons and possible couplings between 
themselves and to nucleons. Following the standard one boson exchange models, we neglect 
meson-meson correlations and treat the interaction between mesons and nucleons perturba- 
tively. The self-interaction of each meson makes the difference to standard models as it is 
taken into account in a closed analytic form and persistently, leading to analytical solutions 



of the nonlinear field equations for each Fourier component of the meson fields in Eq. ( p.ll) . 
For the explicit form and the quantization of these solutions, dubbed as solitary meson 
fields, we refer to our former publication ||^. 

The probability for the propagation of solitary mesons can now be defined as the ampli- 
tude to create an interacting field at some space-time point x which is annihilated into the 
vacuum at y. The momentum space amplitudes of the solitary meson propagator then reads 



k,n 
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where we introduced the dimensionless coupling constants 
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The Feynman propagator 
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uses the mass spectrum 

M^^p = (2pn + l)mp. 

For p = 1/2 one gets the amplitude for scalar fields and p = 1 describes pseudoscalar 
particles. Vector mesons require p = 1 and each term of the sum is multiplied with the 
Minkowski tensor 



-9'" + 



n,v , 



The series ( |2.2|) converges rapidly and in practical calculations it is sufficient to use n < 4. 

The Lorentz invariant normalization D^^l^, which occurs in the propagator ( |2.2| ), is ob- 
tained from the normalization d[^^ of the solitary meson fields by substituting 

i k^". 

2pn + 1 

At this point, we simplify our model. The linear cr-model implies that the nonlinear 
term associated with A2 in Eq. is zero for all mesons despite the scalar a and S 

mesons. Since the former is an effective particle to simulate two-pion exchange and the 
latter contributes little, it is no disadvantage to put A2 = 0, implying = 0, for all mesons 
used. This restriction simplifies our expressions and we notice ai = which allows to drop 
the subscripts from and . 



B. Proper normalization 



The momentum dependent normalization D^f^ of the solitary mesons plays an important 
role and requires a detailed discussion. D^f^ can depend on the four-momentum k'^ and the 
coupling constant The following conditions [0 are imposed: (i) all amplitudes are to 
be Lorentz invariant, (ii) D^^^ is dimensionless, (iii) all Feynman- diagrams are to be finite, 
and (iv) the fields are vanishing for 0. 

The amplitude ( p.2|) has to fulfill on-shell conditions known from renormalization theory 
| T5| . At k"^ = m^, the propagator iP{k'^ ,mp) has to have a pole with residue i. Defining 

(2)/, 2n r_ o 
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these conditions can be met using 
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and 



(2.8) 



From Eq. ( p.2|) it is clear, that iP{k'^ ,mp) readily fulfills Eq. ( p.7|) . The second condition 
demands that D^^^ equals one for /c^ = m^. 

Furthermore, to obtain finite results for all self-energy diagrams involving solitary 
mesons, it is sufficient to choose 

= 0{k% 

for spinless particles, and 

Df = 0{k% 

for vector mesons. In summary, all conditions are met with 



where S denotes the particle spin. With this proper normalization, the solitary meson 
propagator is completely determined and can be applied in a boson exchange potential. 
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C. The scaling law 

In conventional models, meson exchange is described by a product of a Feynman propa- 
gator and an empirical form factor 



(2.10) 



Using the proper normalization ( |2.9| ), the solitary meson propagator was found to resemble 
very closely the expression ( [^.lOD used in the Bonn-B potential 0. This essential result 
permitted us to drop the phenomenological form factors. The astonishing benefit unfolds 
when we make this comparison for all mesons. Doing so, one sees an empirical scaling 
relation for the self-interaction coupling constants 

ap = a^VsTl(—Y . (2.11) 

is the only remaining parameter to describe the full meson dynamics. This reduces 
significantly the number of parameters with respect to Bonn-B potential. 



III. NN PHASE SHIFTS 



In the calculation of A^A^ phase shifts, we use the meson masses of the Bonn-B potential. 
There is some evidence that the ttNN coupling constant should have a value below the 
previously used one g'^/4:TT = 14.4. The first indication came from a Nijmegen analysis |l7 



which suggests f^j^j^ = 0.0745 (which yields gl/4:n = 13.79 with our values for the pion 
and nucleon mass). Additionally, Arndt and co-workers deduced similar value from their 
analysis of tcN scattering [T^. Since we confirmed their result in an independent analysis 
]TP| and intend to apply OSBEP in vrA^ interactions, we fix the ttNN coupling constant to 
the Arndt value 

^ = 13.75. (3.1) 

The parameter and the remaining meson-nucleon coupling constants then yield a total 
number of eight adjustable parameters. 



As in the first analysis 0, we started our fitting procedure with phase shifts |T4 
and deuteron properties, disregarding data. In this case, we are free of Coulomb effects 
and there are more partial waves due to the isoscalar and isovector contributions. In this 
work, the data were also considered. This required to replace (i) the average nucleon mass 
938.926 MeV by the proton mass 938.272 MeV, (ii) the average pion mass 138.03 MeV by the 
7r° mass 134.98 MeV. Additionally, the (Xi coupling constant was reduced by 1.3% from its 
np value. A static point charge Coulomb potential was included using the Vincent-Phatak 
method to calculate the Coulomb distorted hadronic phase shifts |^ 



The final parameter set is listed in Table |. Deuteron properties are very well reproduced 
and are given in Table ||. np phase shifts are shown in Figs.|I| and ^ for single and coupled 
channels, respectively, phases are contained in Fig. ^ We plot the results from the 
Bonn-B, Nijm93, Paris, and OSBEP potential as well as the single energy SM97 analysis. 
All potentials are in close agreement. Differences, with values of several degrees, do exist 
for the ^5*0 phase shifts of which the Paris potential is the worst. This is true over the 



8 



whole energy range. Experimentally, the change of sign for np lies at Tiab = 255.2 MeV ||12|| . 
At 250 MeV, the theoretical values are -2.47 (Paris), -1.72 (Bonn-B), -0.45 (Nijm93), 0.73 
(OSBEP), and 1.03 ±0.84 (SM97). Most striking are the deviations in the P channels. They 
become crucial at energies above 50 MeV and are visible in the observables. There exist a 
large amount of literature about these deviations, but a convincing and final solution has 
not been put forward. In particular, it is known that a potential without vrp correlations 



leads to an overattraction in the P waves |211]. This can be expected to have a major effect 



in pp scattering, since isoscalar D waves are absent. It is surprising that, despite the large 
and consistent database which determines the phase shifts and which is well described by 
the potential models, there are none the less some strong deviations within the model phase 
shifts. 

Our fitting procedure leaves the coupling constants in qualitative agreement with most 
Bonn-B values 0. Differences occur for the ttNN and rjNN coupling constants and the 
tensor to vector ratio n. We use the experimental value 5'^/47r = 13.75, whereas Bonn-B 
uses 14.4 which was preferred in the 1980's. Differential cross section data at backward 
angles support the lower value. SU(3) flavor symmetry relates the irNN and rjNN coupling 
constants to be 

with Of ~0.6-0.65. This yields 0.7 < 5'^/47r < 1.7, consistent with the value in Table |. 
More support for the small value can be found in literature p2|. Increasing the riNN cou- 



pling serves to simulate vrp contributions which are generally absent in one boson exchange 
potentials In the Bonn-B potential, the value g'^/4:7i = 3 is used. 

Another feature of our parameter set is the low tensor to vector ratio k = 3.8 which 
is in close agreement with the vector-dominance value 3.7, to be compared with k = 6.1 
in Bonn-B. This is reconciled by introducing a direct vector coupling of the photon to the 
nucleon pil. We agree with the Nijmegen group that a small nNN coupling constant should 



be aligned with a value of n close to the vector dominance value [T^ . 



IV. OBSERVABLES OF NN SCATTERING 



To obtain observables from phase shifts we follow the notation of Hoshizaki p^. The 



program SAID |12| contains explicitly this option but offers additionally the convention 



of Bystricky et al. [^]. Experimental data with error bars and normalizations together 
with the theoretical phase shifts for Nijm93 and Paris were taken from SAID. Bonn-B we 
calculated ourselves and verified its agreement with published values [0. 

A. np observables 

Altogether, there exist 2719 data points for 13 observables between and 300 MeV. Out 
of 260 possible plots, we selected 19 as representative. They are shown in Figs.0-|^. For each 
measured observable, we plot the theoretical results of OSBEP (full line), Bonn-B (dashed), 
Nijm93 (dotted), and Paris (dash-dotted). Visually, the models are hardly to distinguish. 
This is important in view of quite different phase shifts discussed above and significantly 
different number of adjustable parameters. OSBEP uses about half the parameters of the 
other models. 

However, there are quantitative differences between the models as shown for the 
X^/datum listed in Table|TTl[ The table reflects how the database developed during the 
last years. In the meantime, a number of very precise measurements of differential cross 
sections and polarization observables became available. In particular, the accurate polariza- 
tion data at 183 MeV in Fig. | from the lUCF group [|l^ yield large contributions for the 



Nijm93, Paris, and Bonn-B potential whereas the OSBEP agrees very well with these data. 
Besides that, the differential cross section measurements in Fig. ^ which at large angles are 
sensitive on the vrA^A^-coupling constant, seem to support the low value of 13.75 used in the 
OSBEP potential (full line) rather than the older value of 14.4 which is used in the Bonn-B 
potential (dashed). Therefore, to have a fair comparison, the conventional models should be 
updated to today's database. As yet, the application in np scattering shows that OSBEP is 
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able to describe the data with comparable accuracy as standard A^A^ potentials using eight 
parameters only which lends support for the model of solitary mesons and for the scaling 



law ( |2.11| ) in particular. 

Besides the general excellent agreement in polarization and spin transfer observables 
we stress the high accuracy which is obtained in the description of the np spin-correlation 
parameter A^z at 67.5 MeV, Fig.||, measured by the Basel group In this context, 

Klomp, Stoks, and de Swart |^ argue that a potential which describes A^z at this energy 
does not allow a high ^SDi mixing angle ei at 50 MeV. Their own PWA, including the Basel 
data, yields ei = 2.2° ± 0.5° at this energy. Figure |^ shows that all models considered here 
predict ei slightly below 2°. This must be compared to ei = 2.9° ±0.3°, a value obtained in a 
phase shift analysis based on the Basel data [|25|, and the SM97 value which is 2.53° ± 0.19°. 
We are inclined to follow the arguments of Klomp et al. that these values are too large. A 



similar argument was given by Machleidt and Slaus |27 



The total elastic cross section is very well accounted for by OSBEP over the whole energy 
range, see Fig. 0. The x^/datum equals 9.5 for 319 data points below 300 MeV. This value 
is surprisingly high and not anticipated from Fig. 0. This is mainly due to one dataset only, 
the experiment of Lisowski et al. whose 67 data points are associated with very small 
error bars, which contributes a x^/datum of 39. The remaining 252 atot measurements are 
fitted with x^/datum of 1.6. 



An other quality of fit is obtained by the high-precision A^A^ Bonn-CD [|T0[ and the Rei- 
dlike Nijmegen potentials [|ll|]. They sacrifice the simplicity of the original boson exchange 
potentials and fit each partial wave separately. 



B. pp observables 

Experimental data cover the interval 1-300 MeV. The data below 1 MeV have been dis- 
carded since an assessment of the low-energy data is difficult in the sense that the full 
electromagnetic interaction has to be taken into account, which is very hard to do in a 
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momentum space calculation. Additionally, these data are associated with very small errors 
and it is misleading to include them in a calculation since they can easily distort the 



result After this subtraction, we are left with 1292 data points for 16 observables. In 
Figs. we show 30 plots representative for a total number of 215 possible plots. The 
Paris potential is still a good fit to the data. As the in Table |TT| indicates, the qual- 
ity of OSBEP is not as good as in the case of n-p scattering. This can be traced to the 
overattraction of one boson exchange potentials in P channels which signals the lack of vrp 



correlations in our potential. A more detailed discussion can be found in [^. np data are 
less seriously affected, since isoscalar contributions partly compensate this shortcoming. We 
have noticed that a significant improvement is achieved with an artificially large rjNN cou- 
pling constant, which contradicts the SU(3) flavor symmetry constraint ( |3.2|) . The full Bonn 
potential includes up correlations and neglects i] exchange, putting = 0. The one boson 
exchange approximation Bonn-B simulates the same contributions by using g'^/4:TT = 3. We 
prefer to use a value g^/4:'K = 0.702 which agrees with SU(3) symmetry and rely on a more 
elaborated model, including vrp contributions and A isobars, to provide better P wave phase 
shifts in a future work. 

A pp version of Bonn-B does not exist in literature. We generate a Bonn-B potential 
suitable for pp analysis by substituting the average nucleon and pion mass of the np version 
by the proton and vr" mass, respectively, and including the Coulomb potential into the 
scattering equation. In addition, we refitted the di coupling constant to be gl^/^'n = 8.8235. 
The same prescription was used for OSBEP. Since the main contribution to the large 
comes from the differential cross section in the energy bin 50-150 MeV, we show some of the 
measured cross sections in Fig. ^. It is obvious, that OSBEP and Bonn-B yield almost the 
same results, in some of the figures the two curves cannot be distinguished. We obtain a value 
for the Bonn-B which is larger than Nijm93 and Paris but slightly below OSBEP. This 
is consistent with the enlarged rjNN coupling constant which somewhat compensates the 
overattraction in the P waves. The remaining harm therefore sticks with the approximations 
made concerning the meson-meson correlations whereas the model of solitary bosons and 
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the scaling law find the same confirmation as deduced from np. 



V. OUTLOOK 

With this analysis, we make a comparison of np and pp observables below pion threshold 
with several potential models. The total x^/ datum shows the high standard of all models but 
also some consistent failures. For the one boson exchange potentials, they become obvious 
for P waves and pp differential cross sections above 80 MeV. This shortcoming is well known 
from older analyses but is here confirmed and has its cause in the absence of meson-meson 
correlations. The phenomenological form factors have been consistently replaced by properly 
normalized solitary meson fields which guarantee finite self-energies. An empirical scaling 
law was deduced from comparison with Bonn-B form factors and this rule was used in case 
of np and pp data. This issue permitted the reduction of fit parameters to the meson-nucleon 
coupling constants and one parameter accounting for the meson self-interaction. This study 
serves the purpose to consistently describe all A^A^ data below pion threshold with the claim 
to be highly quantitative but with significantly reduced degrees of freedom in the fits. The 
X^/datum results are listed in Table pTj The OSBEP result is close to the Nijm93 potential, 
whereas the Bonn-B and Paris potential yield considerably larger values. However, both 
models may easily be refitted to improve their with respect to the latest database. 

As our comparison of several potential models and their predictions for observables of 
elastic A^A^ scattering shows, there is little room for improvements or to discern model 



details on-shell. In previous work [Q, we made a strong point that (p, P7) bremsstrahlung, 
triton binding energy, and nucleon-nucleus scattering cannot discern off-shell differences if 
the on-shell amplitudes are equivalent. 

The boson exchange models cannot be extended towards higher energies, the regime of 
meson production, and hadronic excitations, since this requires a genuine QCD dynamics. 
New experimental facilities, such as lUCF, CELSIUS, COSY, and TJNAF provide high- 
quality data and we are seriously considering various potential models suitable for this new 
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domain. Beyond any doubt, this is a subtle problem. Prior to this, it appears interesting to 
investigate the empirical scaling law in more detail and have a look into the boson exchange 
model for t[N scattering. In this context, it is a common problem that the form factor 
parametrization of the A^A^ interaction can not be used in the calculation of nucleon pole 
diagrams ||5|. This may be the reason for the failure of the attempts to gain a consistent 
description of NN , ttN and vrvr interactions. To achieve this goal would lend support for 
the proper normalization of solitary meson fields. The higher order diagrams in the vrA^ 
scattering equations, which need to be regularized by a form factor, are essentially baryon 
self-energy and vertex correction amplitudes. Since the proper normalization was designed 
to yield finite results for these diagrams, it is corollary to work also there. 
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TABLES 



TABLE I. OSBEP parameters. 





IT 


V 


P 


u 




0-1 


6 




Q- 


0- 


1- 


1- 


0+ 


0+ 


0+ 


[MeV] 


138.03^ 


548.8 


769 


782.6 


720 


550 


983 


4tt 


13.75 


0.702 


1.431 


21.07 


14.64 


8.6619^ 


1.259 


= 0.44065 








fp/9p 


= 3.829 





'^Values for the pp potential are 


= 134.9764 MeV and gljA-K 


= 8.5531. 








TABLE IL Deuteron properties 








Bonn-B 


§ OSBEP 


Exp. 


Ref. 


Eb (MeV) 


2.2246 


2.22459 


2.22458900(22) 


m 


fJ-d 


0.8514^ 


0.8524^ 


0.857406(1) 




Qd (fm2) 


0.2783^^ 


0.2698^ 


0.2859(3) 




As (fm-V2) 


0.8860 


0.8805 


0.8802(20) 




D/S 


0.0264 


0.0258 


0.0256(4) 


il 


rRMS (fm) 


1.9688 


1.957 


1.9627(38) 




Pd (%) 


4.99 


4.80 







'^Meson exchange current contributions not included. 



TABLE IIL x^/datum for the OSBEP and several potential models. Data and values for 



the OSBEP, Nijm93 and Paris potential were taken from SAID [12| 



Model 
OSBEP 
Nijm93 
Bonn-B 
Paris 



No. of param. 
~8 
15 
15 
«60 



np 

4.1 

5.6 

12.1 

17.5 



pp 

6.8 

2.2 

5.8'= 

2.3 



Total 

5.0 

4.5 

10.1 

12.6 



^Energy bin 0-300 MeV (2719 data points). 
^Energy bin 1-300 MeV (1292 data points). 
^pp version g'^_^/ATT = 8.8235, see text. 
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FIG. 1. np phase shifts. We show the Arndt SM97 phase shift analysis (circles) compared 
to the potentials Nijm93 (dotted), Bonn-B (dashed), Paris (dash-dotted), and OSBEP (full). 
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20 




21 






FIG. 4. Observables of np scattering. Kinetic laboratory energy is denoted, experimental data 
are taken from |jl2| with notation from We show theoretical predictions from OSBEP (full) 
and Bonn-B (dashed), Nijm93 (dotted), and Paris (dashed-dotted). 
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FIG. 5. Observables of np scattering, notations as in Fig. 
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FIG. 6. Observables of np scattering, notations as in Fig. 
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FIG. 7. Total cross section for elastic np scattering. 
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FIG. 8. Observables of pp scattering. Kinetic laboratory energy is denoted, experimental data 
are taken from with notation from |2^. We show theoretical predictions from OSBEP (full), 
Bonn-B (dashed, see text), Nijm93 (dotted), and Paris (dashed-dotted). 
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FIG. 9. Observables of pp scattering, notations as in Fig. 
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